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ABSTRACT 

Joint  probability  distributions  and  conditional  means 
functions  are  used  to  analyze  turbulent  wind  and  temperature 
data  obtained  over  natural  ocean  waves.   The  data  were 
collected  aboard  R/V  FLIP  during  BOMEX  1969. 

Stable  conditions  prevailed  during  all  time  periods 
considered.   The  emphasis  is  to  use  the  statistical  pro- 
cedures to  describe  the  wave's  influence  on  the  temperature 
fluctuations  at  the  8  meter  level. 

Results  are  compared  with  predictions  given  by  potential 
flow  theory.   The  results  consistently  demonstrate  that  waves 
exert  an  influence  on  the  temperature  fluctuations.   In 
addition,  it  is  observed  that  under  near  neutral  thermal 
stratification  or  prolonged  periods  of  stable  thermal 
stratification  the  temperature  fluctuations  approached  those 
expected  from  potential  flow. 
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I.   INTRODUCTION 

Turbulent  and  wave  induced  temperature  fluctuations  over 
ocean  waves  may  be  described  by  theoretical  models  or  inferred 
from  observations.   These  descriptions  are  necessary  to 
■develop  an  understanding  of  a  wide  class  of  applied  problems 
which  require  a  knowledge  of  sensible  heat  transfer  and  the 
nature  of  temperature  fluctuations. 

Vertical  turbulent  transfers  of  both  sensible  and  latent 
heat  in  the  boundary  layer  are  influenced  hy  modifications 
induced  in  the  overlying  air  mass  by  the  surfaces  below. 
Most  present  theories  and  procedures  for  estimating  these 
transfers  are  based  to  a  large  extent  on  empirical  formulae 
developed  from  data  obtained  over  land.   These  include 
assumed  profiles  for  the  mean  wind,  temperature   and  humidity. 

However,  recent  theoretical  and  observational  studies 
that  will  be  described  in  the  background  section  suggest  that 
transfer  processes  over  ocean  waves  may  be  influenced,  and 
not  necessarily  in  a  simple  manner,  by  the  waves  themselves. 
Therefore,  one  cannot  assume  that  phenomena  which  are  fairly 
well  understood  over  land  are  also  understood  over  water 
until  the  wave's  influences  are  adequately  determined. 

One  applied  problem  in  which  a  better  description  of  the 
temperature  fluctuations  over  waves  is  required  is  the 
optical  (laser)  propagation  problem.   These  studies  require 
a  knowledge  of  the  refractive  index  in  order  to  predict  the 
path  arid  the  scintillation  of  optical  wave  propagata 


Temperature  fluctuation  would  influence  this  index  and 
thereby  optical  wave  propagation.   Optical  wave  propagation 
is  presently  based  on  stability  parameters  which  ignore  the 
influence  of  waves. 

The  purpose  of  this  study  is  to  examine  turbulent  features 
of  the  temperature  structure  over  natural  ocean  waves.   The 
observational  data  were  obtained  aboard  R/V  FLIP  (Floating 
Instrument  Platform)  during  BOMEX  (Barbados  Oceanographic  and 
Meteorological  Experiment)  during  the  period  19  to  28  May 
1969.   The  data  examined  are  from  simultaneous  measurements 
of  velocity  components  (u,v,w),  temperature  (T)  at  fixed 
levels  above  mean  sea  level  and  wave  height  (n) . 

The  analyses  methods  used  in  this  study  are  the  same  as 
those  in  an  investigation  of  velocity  fluctuations  over  waves 
by  Davidson  and  Frank  (1973).   However,  temperature  fluctua- 
tions are  the  parameters  of  interest  in  this  study  whereas 
the  velocity  fluctuations  were  considered  in  the  referenced 
study.   Statistical  procedures  applied  in  these  analyses  are 
joint  probability  density  function  and  conditional  mean 
function.   The  former  will  be  referred  to  as  (JDF)  and  the 
latter  as  (CMF)  in  the  remainder  of  this  thesis.   This  study 
was  the  first  application  of  the  JDF/CMF  analysis  to  tempera- 
ture data  over  ocean  waves.   Results  will  be  discussed  in 
relation  to  observational  studies  by  Volkov  (1969),  Phelps 
and  Pond  (1971)  and  Pond  et  al .  (1971). 

In  addition  to  the  JDF /CMF  statistical  procedures  this 
study  will  include  spectral  results •   Thompson  (1972)    s 
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conducted  a  study  of  momentum  transfer  features  using  the 
same  statistical  procedures  and  the  same  time  periods. 
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XI.   BACKGROUND 

A  study  of  boundary  layer  wave  induced  fluctuations  re- 
quires some  sort  of  physical  theory  connecting  observed 
properties  to  those  we  wish  to  infer.   This  physical  theory 
•may  be  obtained  by  developing  a  model  which  gives  realistic 
results  or  by  observing  real  data  and  describing  what  has 
occurred.   The  results  of  such  studies  will  frequently  allow 
inferences  of  properties  other  than  those  that  are  measurable 
directly  by  instruments.   This  background  discussion  will 
first  consider  some  papers  presenting  theoretical  approaches 
and  then  review  recent  attempts  to  use  real  data. 

A.   THEORETICAL  APPROACH 

With  respect  to  theoretical  predictions  Stewart  (1969) 
has  suggested  that  with  a  stable  thermal  stratification  the 
study  of  turbulent  fluctuations  is  complicated  by  the  fact 
that  motion  due  to  both  turbulence  and  internal  gravity 
waves  is  present  and  the  aspects  we  would  like  to  infer  are 
totally  different  for  each  kind  of  motion.   In  addition  to 
the  existence  of  two  different  types  of  fluctuations,  there 
is  a  nonlinear  coupling  between  waves  and  turbulence. 
Stewart  (1969),  while  looking  for  a  tracer  of  turbulent 
motion,  stated  that: 


it 

■  . 


the  waves  are  capable  of  transporting  momentum 
but  they  are  not  capable  of  transporting  passive 
scalars  per  se." 

Therefore,  he  suggests  using  a  scalar  tracer  transport   n 

order  to  see  turbulent  effects.   This  in  eifect  imp ]     that 
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a  scalar  such  as  temperature  could  be  predicted  by  potential 
flow  theory.   Velocity,  temperature  and  heat  flux  stream- 
lines shown  in  Figure  1  are  those  that  would  occur  under 
stable  temperature  stratification  and  potential  flow.   This 
concept  of  the  coupling  between  a  scalar  temperature  field 
over  a  natural  wave  field  is  in  contrast  to  that  presented 
by  Volkov  (1969)  which  will  be  reviewed  later  in  this  text. 

In  a  theoretical  study,  Witting  (1972)  considered  an 
irrotational ,  plane,  periodic,  linear  progressive  wave  model 
to  study  the  variation  of  surface  temperature  induced  by 
waves  and  thereby  make  it  possible  to  estimate  the  amount  of 
phase  shift  away  from  potential  flow.   He  concluded  that  (1) 
when  the  wave  period  is  sufficiently  small (^  3  sec)  the 
amplitude  of  the  wave  induced  temperature  fluctuation  is 
proportional  to  the  wave  amplitude  and  to  the  square  root  of 
the  wave  period,  (2)  surface  temperature  maxima  lead  wave 
crests  by  one-eighth  of  a  period  when  the  average  surface 
temperature  is  less  than  the  bulk  water  temperature  and  (3) 
surface  temperature  minima  lead  wave  crests  by  one-eighth  of 
a  period  when  the  average  surface  temperature  exceeds  that 
of  the  water  below.   This  theoretical  study  then  recognizes 
wave  induced  fluctuations  in  the  scalar  temperature  field 
and  a  deviation  from  potential  flow. 

B.   OBSERVATIONAL  APPROACH 

Since  theoretical  studies  do  not  agree,  attention  is 
directed  toward  observational  studies.   Phelps  and  Pond  (1971) 
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and  Pond  et  al .  (1971)  examined  data  obtained  aboard  R/V 
FLIP  during  BOMEX  1969.   However,  in  their  study  wave  motion 
was  not  measured  and  no  recognition  is  given  to  wave  induced 
temperature  fluctuations.   According  to  similarity  theory 
predictions,  the  intensities  of  velocity  and  temperature 
fluctuations  should  depend  only  on  the  wind  speed,  the  rough- 
ness of  the  underlying  surface   and  the  thermal  stratifica- 
tion.  In  the  case  of  temperature  fluctuations,  this  depend- 
ence can  be  expressed  as 

2-  =  f(z/D 

where 

T   =   I  .  SI' 
K         uA 

T   =  (T-T)1/2 
u^  -  (-  uw) 

L  =  -  TuA3/Kg  wT' 

K   =  von  Karman ' s  Constant 

W  '   =  W  -  TT 

w   =  observed  vertical  motion 
w   =  mean  vertical  motion 
T'  =  t-T 

T   =  observed  temperature 
T   =  mean  temperature 
Z   =  height  above  mean  sea  level 
Therefore,  the  ratio  aT/TA  should  depend  only  on  Z/L.   Such 
a  dependence  has  been  verified  in  numerous  ov     rid 
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experiments.   In  particular,  the  ratio  has  been  observed  to 
haye  a  value  of  approximately  1  in  near  neutral  conditions. 

Phelps  and  Pond  (1971)  observed  agreement  between  this 
aT/TA  results  and  the  similarity  theory  predictions  with 
respect  to  Z/L.   However,  they  also  observed  that  (JT/T^ 
varied  as  c/uA,  where  c  corresponds  to  the  phase  speed  of 
the  predominant  surface  wave.   Kitaygorodskiy  (1969)  and 
Volkov  (1969)  initially  suggested  that  such  a  relationship 
should  exist.   Because  aT/TA  is  inversely  proportional  to 
the  efficiency  of  the  turbulent  flow  to  transport  sensible 
heat,  the  latter  findings  imply  that  over  a  wave  surface 
there  are  considerable  temperature  fluctuations   which  are 
not  associated  with  heat  transfer.   The  latter,  of  course, 
was  suggested  by  Stewart  (1969)  with  respect  to  wave  motion 
within  a  thermally  stratified  layer. 

The  observed  relationship  between  aT/T^  and  c/uA  appears 
in  Figure  2.   The  curves  mean  square  deviations  are  those 
constructed  by  Volkov  (1969).   This  suggested  curve  will  be 
compared  later  to  results  obtained  in  this  study. 
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Figure  2.   aT/ | T ^ |  vs  c/u A  According  to  Monin  Obukhov 
Similarity  Theory. 


Phelps  and  Pond  (1971)  used  data  similar  to  Pond  et  al. 
(1971)  to  describe  the  similarities  and  differences  between 
temperature,  humidity,  and  velocity  fluctuations.   They  found 
that  over  a  wide  frequency  range  the  temperature  waveform 
leads  the  vertical  velocity  and  humidity  waveform   by  20-40 
degrees,  but  they  attributed  this  to  instrument  effects  and 
R/V  FLIP's  motion.   Again  it  should  be  noted  that  no  wave 
data  were  available  for  their  experiment.   They  also  con- 
cluded that  humidity  and  vertical  velocity  are  basically  in 
phase.   No  humidity  data  were  available  for  the  present  study 

In  another  observational  study  Volkov  (1969)  reduced 
western  Mediterranean  data  according  to  similarity  theor" 
In  order  to  compare  characteristics  of  turbulence  in 
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atmospheric  boundary  layer  above  an  agitated  sea  surface  with 
measurements  in  the  boundary  layer  above  dry  land.   He  states 
that  the  main  difficulty  in  interpreting  atmospheric  measure- 
ments above  the  sea  is  the  presence  of  a  mobile  fluid  air- 
water  interface  which  has  an  irregular  variation  as  a  result 
of  wind  action  and  airflow-wave  feedback  processes.   The  feed- 
back is  due  to  the  surface  waves  disturbing  the  airflow 
around  them  and  forming  waves  in  the  air  which  interact  with 
the  velocity  gradient  of  the  basic  flow.   The  wave  induced 
motion  in  the  air  flow  causes  some  energy  to  be  transferred 
from  the  basic  flow  to  the  waves.   These  processes  also 
effect  the  fluctuation  of  a  scalar  such  as  temperature  or 
moisture.   In  particular  Volkov  (1969)  suggests  that  an 
observed  temperature  (T)  above  a  wave  surface  should  be  made 
up  of  the  mean  temperature  (T(z.)),  turbulent  fluctuations 
(T * (x,y , z , t ) )  and  wave  induced  fluctuations  (T"(x,z)). 
Therefore : 

T  =  T(z)  =  T"(x,z)+  T' (x,y,z,t) 

Therefore,  of  special  significance  in  studying  the 
structure  of  turbulence  in  the  atmospheric  boundary  layer 
above  the  agitated  sea  surface  are  measurements  of  the  fluc- 
tuations of  velocity  (u,w),  temperature  (T)  and  emphatically, 
wave  height  (n).   It  should  be  noted  that  Phelps  and  Pond 
(1.971)  and  Pond  et  al  .  (1971)  did  not  have  simultaneous 
measurements  of  wave  height  included  in  their  data. 

By  the  use  of  spectral  analysis  Volkov  (1969)  was  al  Le 
to  conclude  that  the  spectra  of  fluctual       al      tl 
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agitated  sea  surface  are  due  to  two  types.   One  type  of 
fluctuation  corresponds  to  those  smaller  than  the  charac- 
teristic scales  of  sea  swell  and  is  close  to  the  correspond- 
ing interval  of  the  turbulence  spectrum  for  the  air  layer 
above  land.   The  other  type  corresponds  to  fluctuations  of 
the  order  of  the  characteristic  size  of  the  surface  waves 
and  is  typlif  ied  by  the  superposition  of  wave  disturbances  on 
the  random  turbulent  field  of  the  fluctuations.   The  latter 
type  was  shown  in  Figure  1  as  potential  flow. 

Davidson  and  Safley  (1973)  in  another  observational 
study  described  the  occurrence  of  micro-thermals ,  large 
positive  temperature  fluctuations,  and  provided  arguments 
supporting  the  concept  of  wave  influence  on  micro-thermals. 
A  review  of  the  literature  associated  with  wind-wave  coupling 
was  provided  by  Davidson  and  Safley  (1972)  who  strongly 
suggest  that  the  wave's  influence  could  contribute  to  micro- 
thermal  formation  in  the  case  of  two  proposed  mechanisms. 
These  mechanisms  were  the  existence  of  wave-related  motion 
in  airflow  which  is  too  organized  to  be  viewed  as  turbulent 
and  instantaneous  distortions  of  the  wind  shear  over  differ- 
ent phases  of  the  surface  wave.   Velocity  spectral  results 
from  a  Lake  Michigan  study  were  used  to  show  that  the  waves 
could  have  played  a  role  in  the  observed  micro-thermal 
phenomena  due  to  wave  related  vortices  in  the  mean  shear. 
However,  large  thermal  instability  associated  with  the 
periods  used  in  the  referenced  study  did  not  allow  a  detail 
description  of  the  precise  role  of  the  waves. 
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On  the  basis  of-  the  results  of  Witting  (1972),  Volkov 
(1969)  and  Davidson  and-Safley  (1972)  just  described,  there 
is  reason  to  expect  that  the  wave  surface  has  a  considerable 
influence  on  all  of  the  statistical  characteristics  of  tur- 
bulence in  the  near  surface  layer.   Volkov  (1969)  states  this 
fact  in  the  following  manner; 

"We  therefore  succeed  in  finding  some  typical  features 
and  regularities  of  this  influence  by  the  complex 
method  of  investigating  the  processes  of  interaction 
between  the  atmospheric  boundary  layer  and  the  agitated 
sea  surface .  " 

The  present  study  is  an  attempt  to  describe  the  above 

processes  through  the  use  of  actual  data  and  statistical 

procedures  to  test  concepts  put  forth  by  Volkov  (1969)  whose 

results  were  obtained  by  semi-empirical  arguments. 
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III.       DESCRIPTION  AND  MEASUREMENT  OF  DATA 

A.  LOCATION  OF  R/V  FLIP- 

The  data  were  from  measurements  taken  by  University  of 
Michigan  personnel  during  the  BOMEX  experiment  of  1969.   The 
data  consisted  of  velocity,  temperature,  vertical  motion 
and  wave  height  measurements  at  three  and  eight  meter  levels. 
The  data  examined  were  from  19,  26,  and  28  May  1969.   During 
this  period  R/V  FLIP  was  located  northeast  of  Barbados, 
Figure  3A. 

FLIP's  position  into  the  wind  was  held  by  a  tug  attached 
with  a  bridle  and  approximately  800  meters  of  line.  The  tug 
then  steamed  slowly  downwind. 

B.  SENSOR  MOUNTING  AND  INSTRUMENTATION 

Instrument  location  relative  to  R/V  FLIP  is  shown  in 
Figure  3B.   As  indicated,  all  sensors  were  mounted  on  the 
vertical  mast . 

Temperature  fluctuations  were  measured  by  a  resistance 
thermometer,  Figure  4A.   The  sensor  consisted  of  a  30  ohm, 
0.00015  inch  diameter  tungsten  wire,  0.16  inches  long  as 
part  of  a  constant-current  circuit  with  a  small  filament 
current  (about  2ma)  to  serve  as  a  resistance  thermometer. 
The  frequency  response  of  the  system  is  determined  primarily 
by  the  thermal  inertia. 
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(A) 


(from  Portman  et  al.  1970) 


(B) 


-%m 


-Zm 


Figure  3.   Arrangements  during  BOMEX  (A)  Geo- 
graphic position,  (B)  Sensor  location  on  board 
K/V  FLIP. 
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A  block  diagram  of  the  bridge  circuit  is  shown  in  Figure 
4B.    An  unbalance  of  the  bridge  is  caused  by  temperature 
and  thereby  a  sensor  resistance  change. 

System  noise  and  sensitivity  restricted  spectral  descrip- 
tion of  temperature  to  frequencies  below  20  Hz .   Since 
sensible  heat  flux  is  associated  with  frequencies  below  10  Hz 
(Deacon  1959),  the  response  was  sufficient  for  this  aspect 
of  the  interpretations.   Furthermore,  of  particular  interest 
is  the  frequency  response  near  the  frequency  of  the  dominate 
surface  wave,  approximately  0.1  Hz. 

-  Velocity  measurements  were  made  with  a  cons tant-tempera- 

2 
ture  hot-film  anemometer   system  using  three  sensor  probes. 

A  discussion  of  this  equipment  is  provided  by  Portman  et  al . 
(1970)  . 

Visual  estimates  of  wave  height  and  10  meter  wind  speed 
are  shown  in  Figure  5.   All  periods  considered  are  charac- 
terized by  stable  thermal  stratification  with  swell  type 
waves  propagating  through  the  area. 


Flow  Corporation  System,  Series  900 

2 

Thermo-Sys tem  Incorporation,  Model  1054B 


23 


(A) 


BRIDGE- 
AMPLIFIER] 


MONITOR 


SUPPRESSOR -AMPLIFIER 


-/VW~ 


■A/W"*-" 


W— i 


(B)  (from  Davidson  1970) 

Figure  4.   Temperature  system  (A)  Diagram  of  the 
temperature  probe,  (B)  Bridge-amplifier  system 
(constant  current). 
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IV.   STATISTICAL  CONSIDERATIONS  AND  PROCEDURES 

A.   SUMMARY  OF  JDF/CMF  PROCEDURES 

Joint  probability  density  function  (JDF)  and  conditional 
mean  function  (CMF)  statistical  procedures  used  in  this  study 
are  based  on  descriptions  by  Holland  (1968,  1972).   The  same 
procedures  were  applied  by  Frank  (1971) t    Davidson  and  Frank 
(1973)  and  Thompson  (1972)  to  small  scale  fluctuations  over 
waves . 

A  JDF  allows  simultaneously  observed  pairs  of  values  to 
be  represented  by  contours  on  a  two  dimensional  array.   A 
trivariate  statistical  relationship  is  obtained  by  computing 
a  conditional  mean  of  a  third  variable  as  a  function,  CMF, 
of  the  two  independent  variables  contributing  to  the  JDF. 
These  statistics  are  determined  for  variables  that  are  nor- 
malized by  their  standard  deviation. 

A  JDF/CMF  array  appears  in  Figure  6A  in  which  both  the 
ordinate  and  abscissa  are  given  in  standard  deviation  (a) 
units.   The  array  is  centered  on  the  mean  value,,  which  is 

zero,  of  the  two  JDF  variables.   The  bottom  number  on  each 

3 
cell  represents  the  JDF  (the  probability  x  10   of  the  joint 

occurrence  of  the  pair) .   A  nine  point  smoother  was  used  to 
remove  the  perturbations  viewed  as  being  physically  and 
statistically  insignificant.   Statistical  procedures  neces- 
sary to  obtain  final  JDF/CMF  results  are  shown  in  Figure  7 
and  will  be  described  briefly  in  the  following  paragraph.  . 
A  complete  discussion  of  the  theory  associated  with  t 
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Figure  6.   JDF/CMF  Array  (A)  Example  of  printed  JDF/ 
CMF  array,  (B)  Octant  and  wave  amplitude  class 
boundary  used  for  phase  amplitude  determination. 
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NORMALIZATION  AND  COMPUTATION 
Of    STATISTICS  FOR  THE  VARIABLES 

BANDPASS 
FILTERING 


COMPUTATION  OF  JDF/CMF, 
APPLICATION  OF  A  9  POINT  SMOOTHER, 
COMPUTATION  OF  STATISTICAL  EDDY  PERIODS, 

AND  PRINT  OUT  OF  ARRAY  3 


COMPUTATION  OF  PHASE 
AMPLITUDE  INFORMATION  USING 
PHASE  ANGLE 


CLASSIFY  INTO  LARGE  AND 
SMALL  AMPLITUDE  WAVES 


PLOTTING  OF  JDF/CMF 

VARIABLES  AND  PHASE 

AMPLITUDE 


Figure  7.   Block  diagram  of  computer  steps  to 
produce  a  smooth  JDF/CMF  array. 
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procedures  was  presented  by  Holland  (1968,  1972),  and  a  more 
complete  discussion  of  its  application  in  this  study  is  pro- 
vided by  frank  (.1971)  . 

B.   BANDPASS  FILTER 

The  bandpass  filter,  shown  as  step  1,  was  an  inverse 
transform  filter.   This  was  applied  to  isolate  the  signifi- 
cant frequency  band  (0.1  Hz  to  0.8  Hz)  associated  with  the 
ocean  waves.   Its  response  curve  Is  shown  in  Figure  8. 


1.0- 
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/   X 
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I        I 
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n  (HZ)   ' 
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(from  Frank  1971) 


Figure  8.   Response  curve  for  inverse  transform  band- 
pass filter . 


C.   NORMALIZATION  OF  DATA  AND  COMPUTATION  OF  STATISTICS 

Normalization  of  data,  shown  as  step  2  in  the  flow  di- 
gram, was  performed  as  described  in  subsection  A. 
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D.  JDF/CMF  COMPUTATION 

The  JDF/CMF  array,  shown  as  step  3  in  Figure  7,  is  de- 
scribed in  the  beginning,  of  this  section.   Procedures  for 
this  were  described  by  Frank  (.1971)  and  Thompson  (1972). 

E.  STATISTICAL  EDDY  PERIOD  COMPUTATION 

The  "statistical  eddy  period"  computed  in  step  3,  was 
suggested  by  Holland  (1968,  1972)  as  an  approach  to  reveal 
if  the  amplitude  of  the  dependent  variable  such  as  tempera- 
ture varies  nonlinearly  with  the  amplitude  of  the  independent 
variable  such  as  wave  height.   Holland  defined  a  "statistical 

eddy  period"   T   =  2ir  — V 
J  o       ay 

for  any  fluctuating  meteorological  variable  y.   ay  is  the 
standard  deviation  of  the  meteorological  variable  and  ay  is 
its  derivative. 

F.  PHASE  AMPLITUDE  COMPUTATION 

Phase  amplitude  computation,  shown  as  step  4  in  the 
flow  diagram,  was  determined  from  a  JDF  coordinate  system 
with  (n)  as  ordinate  and  its  derivative  (r))  as  the  abscissa. 
This  allowed  presentation  of  the  phase  of  the  variable  with 
respect  to  the  wave.   The  phase  angle  was  defined  by  an 
angle  measured  (by  octants)  counterclockwise  (CCW)  from  the 
horizontal  axes  in  the  0*1,  n)  coordinate  system.   This  is 
shown  in  Figure  6B.   A  complete  discussion  of  this  may  be 
found  in  Holland  (1968,  1972). 

The  phase  amplitude  procedures  outlined  in  the  above 
references  were  performed  with  respect  to  the  waves  foi  che 
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Figure  9.   Example  of  phase  amplitude  results  for 
potential  flow  case. 
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variables  u,  w,  T  and  wT .   Prior  to  summarizing  by  octants 
and  plotting,  the  results  were  separated  into  wave  amplitude 
classes  defined  by  regions  on  the  JDF  array.   Three  classes 
were  considered,  those  being  small  amplitude  (less  than  1.5 
C)  ,  large  amplitude  (greater  than  1.5<j)  and  the  sum  of  the 
first  two.   This  yielded,  therefore,  phase  relationships 
with  respect  to  small,  large  and  all  surface  waves. 

Summarized  phase  amplitude  results  for  each  octant,  as 
described  in  Figure  6B,  are  shown  by  horizontal  bars  that 
appear  in  Figure  9.   Smooth  lines  were  drawn  through  each 
set  of  horizontal  bars  such  that  areas  above  and  below  the 
resulting  curves  were  equal. 

G.   SUMMARY  OF  ANALYSIS  PROCEDURES 

The  previous  discussion  of  JDF/CMF  analysis  procedures 
indicate  that  joint  probability  results  are  well  suited  for 
identification  of  fluctuations  associated  with  waves.   This 
type  of  result  is  particularly  useful  because  it  will  reveal 
non-linear  properties  in  fluctuations.   These  properties 
were  alluded  to  by  Volkov  (1969)  and  Davidson  and  Frank 
(1973)  . 

As  will  be  shown  in  the  presentation  of  results,  some  of 
the  most  significant  interpretations  will  be  with  respect  to 
the  phase  amplitude  relations.   These  results  do  not  repre- 
sent phase  amplitude  information  as  revealed  by  simple  linear 
expansion  of  the  fluctuations,  i.e.  spectral  analysis,  but 
rather  the  actual  nonlinear  average  fluctuations  over  d  f f er- 
ent  portions  of  the  average  wave. 
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V.   GENERAL  WIND,  WAVE,  AND  THERMAL  CONDITIONS 

Data  from  four  separate  measurement  periods  were  con- 
sidered.  These  periods  were  19  May  and  27  May  and  two  on 
28  May  1969.   They  will  be  referred  to  as  periods  1  through  4 
respectively  in  the  following  discussion.   The  times  for  each 
period  appear  in  Table  1  and  Figure  5A  along  with  information 
on  the  mean  ten  meter  wind,  the  significant  wave  height,  the 
heat  flux  (H)   and  Z/L. 

Wave  histories  in  Figure  5A  indicate  that  equilibrium 
conditions  prevail  with  respect  to  the  waves  during  all 
periods.   Sea  conditions  were  generally  dominated  by  swell. 
Period  2  had  the  lowest  significant  wave  height,  1.3  meters, 
while  the  other  three  periods  had  significant  wave  heights 
of  about  2.6  meters.   The  wind  speed  at  the  8  meter  level 
ranged  from  9.5  meters / second  for  period  A  to  11.8  meters/ 
second  for  period  2. 

The  existence  of  equilibrium  conditions  with  respect  to 
atmospheric  and  sea  conditions  for  all  periods  is  evident, 
perhaps,  by  the  fact  that  the  relationship  between  cT/TA  and 
c/u^  are  comparable  to  those  observed  by  Volkov  (1969). 
aT/T^  vs  c/u^  results  from  those  data  are  along  the  horizontal 
section  of  Volkov's  curve,  Figure  10.   Volkov  (1969)  has 
suggested  that  this  range  of  values  corresponds  to  conditions 
with  surge  or  nondeveloping  sea.   The  curve  shown  in  Figure 
10  is  a  replot  of  Figure  2. 
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VI.   PRESENTATION  OF  RESULTS 

A.   SPECTRAL  AND  "STATISTICAL  EDDY  PERIOD"  RESULTS 

Temperature  spectra  for  periods  1,  3  and  4  appear  in 
the  log-log  format  of  Figure  11.   The  spectral  slope  on 
these  plots  at  high  frequencies  approaches  a  -5/3  value  for 
all  three  periods.   The  -5/3  slope  is  expected  for  the 
inertial  subrange  of  a  passive  scalar. 

Pairs  of  temperature,  wT ,  and  wave  spectra  for  periods 
1,  3  and  4  are  presented  in  linear-log  profiles  in  Figure  12, 
13  and  14.   The  left  side  of  each  Figure  is  the  spectrum  at 
the  beginning  of  each  period  while  the  right  side  is  the 
spectrum  occurring  near  the  end. 

The  most  significant  characteristic  of  these  spectra  is 
the  absence  of  any  noticeable  peak  in  the  temperature  spectra 
at  the  predominant  frequency  of  the  wave.   Therefore,  from 
an  examination  of  spectra  alone  it  might  be  mistakenly 
stated  that  ocean  waves  do  not  influence  a  passive  scalar 
such  as  temperature.   However,  statistical  results  will  be 
presented  to  show  that  temperature  fluctuations  are  indeed 
influenced  by  wave  motion  below. 

Spectral  results  as  presented  above  were  not  available 
for  period  2 . 

The  "statistical  eddy  period"  was  computed  as  described 
in  section  4  and  the  results  are  shown  in  Table  2.   The 
statistical  eddy  period  for  the  temperature  waveform  is 
always  less  than  that  of  the  ocean  wave.   This  stror.; 
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suggests  that  potential  flow  theory  will  not  accurately 
predict  the  air  sea  coupling  in  the  boundary  layer. 

B.   JDF/CMF  PRESENTATION  FOR  POTENTIAL  FLOW  CASE 

The  principal  comparisons  in  the  discussion  of  results 
will  be  with  respect  to  potential  flow  theory.   These 
features  are  based  on  simple  vertical  displacement  of  a 
thermally  stratified  air  column  by  the  underlying  progressive 
waves.   If  potential  flow  theory  properly  describes  the 
coupling  of  the  air  column  to  the  waves,  then  fluctuations 
in  temperature  (T)  ,  velocity  (u)  ,  vertical  motion  (w)   and 
heat  flux  (wT)  observed  by  a  fixed  sensor  would  appear  as 
shown  in  Table  1.   JDF/CMF  results  would  be  those  appearing 
in  Figure  15  and  phase  amplitude  results  would  be  those 
appearing  in  Figure  9. 

A  general  description  of  the  JDF  patterns  for  potential 
flow  case  follows : 

1.  Uniform  circles;  90°  phase  shift  between  JDF  variables 

2.  Ellipse;  orientated  135°-315°,  JDF  variables  are  in 
phase 

3.  Ellipse;  orientated  45°-215°,  180°  phase  shift  in 
JDF  variables. 

Dashed  lines  in  all  JDF/CMF  results  correspond  to  the 
zero  isoline  of  the  CMF.   The  orientation  of  various  CMF 
patterns  in  each  panel,  as  shown  in  Figure  15,  is  examined 
closely  in  the  following  interpretations. 
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Log  (f)  Hz 


Figure  11.   Temperature  spectra  for  periods  1,  3  and  4 
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Log  (f)  Hz 

Figure  12.   Temperature,  wT ,  and  wave  spectra  at  (a)  be- 
ginning and  (b)  end  of  period  1. 
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Log  (f)  Hz 

Figure  13.   Temperature,  wT,  and  wave  spectra  at 
(a)  beginning  and  (b)  end  of  period  3, 
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Figure  14.   Temperature,  wT ,  and  wave  spectra  at 
(a)  beginning  and  (b)  end  of  period  4. 
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1.  Panel  A 

Temperature  and  vertical  motion  are  shown  as  the  JDF 
and  wave  height  as  CMF  in  this  panel.   When  temperature  is 
maximum  positive  and  wave  height  is  maximum  negative.   Like- 
wise, the  w  maxima  occur  at  every  other  inflection  point  of 
the  wave . 

2.  Panel  B 

In  this  panel  the  CMF  zero  isolines  separates  the 
positive  and  negative  phase  of  the  temperature  waveform. 
The  JDF  is  the  wave  height  and  its  derivative.   The  octants 
to  be  referred  to  were  defined  on  page  30.   This  figure  may 
be  interpreted  in  a  CCW  direction  starting  in  octant  1. 
Temperature  may  then  be  interpreted  to  show  that  when  the 
wave  height  is  negative  the  temperature  is  positive  and  vice 
versa.   To  state  it  another  way,  this  panel  represents  a 
180°  phase  shift  between  the  ocean  wave  and  temperature  wave 
forms . 

3 .  Panel  C 

The  CMF  zero  isoline  in  this  panel  separates  positive 
and  negative  vertical  motion.   The  JDF  is  again  determined 
from  the  wave  height  and  its  derivative.   Interpretations  of 
this  panel  in  a  CCW  direction  from  octant  1  indicates  that 
if  the  wave  height  is  positive,  w  goes  from  positive  to  nega- 
tive.  This  indicates  a  90°  phase  shift  between  ocean  waves 
and  w  wave  forms  during  potential  flow. 
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4  .   Panel  p 

The  u  and  T  are  presented  in  the  JDF  and  the  CMF  is 
for  wave  height.   In  pctant  1  and  2  u  and  T  are  positive 
over  the  waves  minima  while  in  octant  5  and  6  u  and  T  are 
negative  over  the  wave  maxima. 

5  .   Panel  E 

Heat  flux,  wT ,  is  the  CMF  and  wave  height  and  its 
derivative  are  the  JDF  in  this  panel.   Beginning  in  octant  1, 
interpretation  of  this  figure  will  describe  heat  flux  in 
potential  flow.   During  the  positive  phase  of  the  wave  wT 
completes  a  full  negative  and  positive  cycle  with  similar 
change  for  the  negative  phase  of  the  wave.   This  indicates 
a  45°  phase  shift  with  respect  to  the  wave. 

6  .   Panel  F 

The  CMF  in  this  panel  is  u  and  the  JDF  is  wave 
height  and  its  derivative.   Beginning  in  octant  1,  interpre- 
tation of  the  panel  shows  that  when  the  phase  of  the  wave  is 
positive  the  u  waveform  is  negative  and,  therefore,  a  180° 
phase  shift.   This  is  the  same  phase  relationship  shown  for 
temperature  in  panel  B.   Attention  will  now  be  turned  to  ard 
interpretation  of  results  obtained  from  the  four  periods 
previously  described. 

C.   PERIOD  1  RESULTS 

Panels  A  and  B  of  the  JDF/CMF  results,  Figures  16,  17, 
18  and  19,  will  be  considered  in  detail  for  the  following 
discussion.   A  cursory  examination  of  the  JDF/CMF  patter  s 
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in  panel  A  and  B  of  the  above  figures,  in  comparison  to 
those  in  Figure  15  revealed  that  processes  associated  with 
potential  flow  are  not  the  only  ones  that  occurred. 
Specifically  in  Figure  16A  there  is  a  135°-315°  orientation 
of  the  elliptical  JDF  pattern  and  the  axis  of  the  ellipse 
appears  curvilinear.   The  elliptical  pattern  described  above 
occurred  in  all  remaining  subsets  of  this  period,  Figure  17A, 
18A  and  19A,  and  during  all  periods  considered  in  this  study. 
This  indicates  approximately  a  180°  phase  shift  in  w  and  T 
instead  of  the  90°  observed  in  Figure  15  when  potential  flow 
was  considered.   The  90°  CCW  rotation  of  the  CMF  wave  height 
pattern  indicates  that  in  general  a  w  minima  and  T  maxima 
occur  over  the  wave  crest. 

The  latter  phase  relationships  are  also  apparent  in  the 
phase  amplitude  results  in  Figure  20.   For  potential  flow 

■ 

temperature  maxima  occurred  over  the  wave  trough  and  w  was 
90°  out  of  phase  with  respect  to  the  wave.   A  number  of 
smaller  perturbations  appeared  in  the  T  waveform  of  Figure 
20.   They  could  be  due  to  the  turbulent  contribution  sug- 
gested by  Volkov  (1969)  . 

A  possible  physical  interpretation  of  the  phase  relation- 
ship noted  above  could  be  that  at  a  fixed  observation  level 
the  waves  are  causing  a  transfer  of  warmer  air  downward  and 
colder  air  upward  thereby  inducing  a  net  heat  transfer  down- 
ward.  This  was  not  apparent,  however,  in  the  spectral 
results,  Figure  12,  13  and  14,  where  it  was  observed  that 
significant  spectral  peaks  do  not  appear  in  T  or  wT  at    e 
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predominate  frequency. of  the  wave  spectrum.   Temperature 
and  w  JDF  and  r\    CMF  relationships  throughout  the  remainder 
of  Period  1  appear  in  panel  A  of  Figures  17,  18  and  19. 
The  CMF  patterns  in  panel  A  of  the  above  figures  are  no 
longer  rotated  90°  CCW  from  potential  flow  but  are  approxi- 
mately 45°.   The  45°  rotation  is  most  evident  in  the  large 
wave  class  region,  greater  than  1.5a,  of  the  JDF  array. 
Figure  18A  is  an  excellent  example  of  this  change  in  rotation. 
An  evaluation  of  the  JDF/CMF  results  of  Figure  18  in  conjunc- 
tion with  phase  amplitude  results  in  Figure  22  reveals  that 
W  and  T  are  approximately  45°  out  of  phase.   Potential  flow 

results  had  indicated  a  90°  phase  shift.   The  temperature 

i 
maxima,  considering  only  the  larger  crests,  are  almost  in 

phase  with  the  wave  trough  which  indicates  a  tendency  toward 
potential  flow. 

Changes  in  thermal  stratification  throughout  the  period 
may  have  been  a  factor  in  these  results.   The  thermal  strati- 
fication changed  from  near  neutral  to  stable  and  returned  to 
near  neutral  during  the  course  of  the  period.   The  Z/L 
values  changed  from  -0.03  to  -0.20  and  ended  at  -0.08.   There- 
fore, with  respect  to  stability  and  the  above  results  one  can 
say  that  under  greater  stability  there  was  a  tendency  toward 
potential  flow  as  shown  in  Figure  18A  and  22,  and  there  were 
fewer  perturbations  on  the  basically  sinusoidal  temperature 
waveform  as  shown  in  Figure  22.   Furthermore,  it  was  observed 
in  Figure  19A  and  23  that  as  thermal  stratification  returned 
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toward  neutral  the  JDF/CMF  results  deviate  more  from  poten- 
tial flow  conditions  than  was  observed  in  Figure  18A  and  22 
discussed  above. 

Another  feature  apparent  in  panel  A  of  Figures  16,  17 
and  18  is  that  the  large  wave  class  results  near  the  outside 
.of  the  array  return  to  potential  flow  predictions  more  so 
than  the  small  wave  class  results. 

Panel  B  of  Figures  16,  17,  18  and  19  provided  information 
on  the  phase  relation  of  temperature  with  respect  to  the  wave 
Corresponding  phase  amplitude  summaries  for  these  JDF/CMF 
arrays  appear  in  Figures  20,  21,  22  and  23.   JDF/CMF  results 
in  Panel  B  appear  to  be  more  sensitive  to  changes  in  thermal 
stratification  than  those  in  panel  A.   For  example,  an  exam- 
ination of  Figure  16B  reveals  positive  temperature  over  the 
wave  crests  while  a  similar  comparison  of  Figure  17B  reveals 
negative  temperature  over  the  wave  troughs.   Similarly  an 
examination  of  Figure  18B  reveals  negative  temperature  over 
the  wave  crests  as  the  period  becomes  more  stable.   In 
general  the  JDF/CMF  and  phase  amplitude  patterns  became  less 
well  defined  as  the  period  progressed  and  thermal  stratifi- 
cation became  more  neutral.   As  an  example,  compare  Figure 
18B  with  19B. 

Again  the  small  wave  class  results  differ  more  from 
potential  flow  prediction  than  the  large  wave  class  results. 

Vertical  action  in  association  with  the  ocean  waves  also 
appears  in  panel  C  of  Figures  16,  17,  18  and  19.   This  series 
of  panels  cannot  be  readily  interpreted  for  phase  rela   ms 
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between  the  w  waveform  and  the  ocean  waves.   Summarized  phase 
amplitude  results  for  w  and  the  waves  appear  in  Figures  20, 
21,  22  and  23.   The  results  of  these  figures  also  show  that 
there  is  a  w  maxima  ahead  of  the  wave  crest.   For  the  poten- 
tial flow  case,  there  would  be  a  90°  lead. 

JDF  results  for  the  joint  occurrence  of  u  and  T  with 
respect  to  the  CMF  of  the  wave  height  appear  in  panel  D  of 
Figures  16,  17,  18  and  19.   The  045°-215°  orientation  of  the 
JDF  ellipse  indicates  u  and  T  are  positively  correlated  as 
predicted  by  potential  flow.   However,  the  curvilinear  shape 
of  the  axis  of  the  ellipse  represents  a  deviation  from 
potential  flow  theory.   The  CMF  zero  isoline  for  wave  height 
in  panel  D  is  rotated  45°  CW  from  the  potential  flow  orienta- 
tion.  This  indicates  a  phase  shift  of  both  u  and  T  with 
respect  to  the  wave.   Summarized  phase  amplitude  results  for 
u  and  T  appear  in  Figures  20,  21,  22  and  23. 

D.   PERIOD  2  RESULTS 

Thermal  stratification  was  slightly  stable  at  the  begin- 
ning of  the  period  and  became  near  neutral  as  the  period 
progressed  as  indicated  by  Z/L  values  changing  from  -0.05  to 
-0.01.   Only  the  first  two  subsets  of  the  form  will  be  dis- 
cussed for  this  period  because  all  subsets,  in  contrast  to 
those  just  described  for  period  1,  are  quite  similar.   Sub- 
set 3  and  4  results  appear  in  Appendix  A  in  Figures  32,  33, 
34  and  35. 

In  general  results  reveal  deviations  from  potential   low 
conditions.   This  can  be  readily  noted  by  comparing  patterns 
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Figure  16.   JDF/CMF  Results  for  period  1  subset  1 
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Figure  17.   JDF/CMF  Results  for  period  1  subset  2 
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Figure  18.   JDF/CMF  Results  for  period  1  subset  3 
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Figure  19.   JDF/CMF  Results  for  period  1  subset  4 
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Figure  20.   Phase  amplitude  results  for  period  1 
subset  1(a)  small  waves,  (b)  large  waves,  and 
(c)  all  waves. 
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Figure  21.   Phase  amplitude  results  for  period  1  subset  2 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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Figure  22.   Phase  amplitude  results  for  period  1  subset  3 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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Figure  23.   Phase  amplitude  results  for  period  1  subset  4 
(a)  small  waves,  (b)  large  waves 


and  (c)  all  waves 
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in  Figures  24  and  25.  with  results  in  Figure  15.   Specifi- 
cally, the  CMF  zero  isol.ine  in  Figure  24A  has  a  45   CCW  rota- 
tion from  potential  flow  orientation,  Figure  15A.   A  similar 
rotation  was  observed  in  Figure  18A  which  corresponds  to  the 
more  stable  conditions  in  Period  1. 

The  T  and  w/JDF  results  are  elliptical  in  Figure  24A  as 
in  Period  1  and  has  a  315  -135   orientation.   This  again 
suggests  a  180   phase  shift  between  T  and  w  in  contrast  to 
the  90   predicted  by  potential  flow  theory  shown  in  Figure 
15. 

Relationships  between  features  in  this  period  and  the 
thermal  stratification  can  be  further  examined  by  a  broad 
comparison  of  similar  JDF/CMF  panels  in  Figure  24  and  25. 
The  patterns  in  Figure  24  appear  to  have  fewer  small  scale 
perturbations  than  those  in  Figure  25.   Phase  results  for 
this  period  are  summarized  in  the  phase  amplitude  results  in 
Figures  26  and  27.   It  can  be  noted  that  the  results  in  Fig- 
ure 27  in  contrast  to  the  JDF/CMF  results  in  Figures  25  are 
more  readily  examined  even  though  small  scale  perturbations 
occur.   An  example  is  the  temperature  waveform  in  Figure  27. 
Thermal  stratification  was  near  neutral  during  the  second 
subset  corresponding  to  Figure  27.   This  deviation  from  a 
sinusoidal  waveform,  as  for  Period  1,  is  considered  to  be 
due  to  turbulence. 

With  respect  to  different  wave  classes,  small  wave 
classes  appear  to  approach  potential  flow  features  more  than 
the  larger  wave  classes.   This  feature  is  shown  by  th     ve 
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height  CMF  zero  isoline  in  Figure  24A.  This  is  in  contrast 
to  their  tendency  to  deviate  from  potential  flow  as  was  ob- 
served in  Period  1  during  more  thermally  stable  conditions. 

A  feature  that  did  not  appear  in  Period  1  but  appears  in 
this  period  is  the  tendency  for  the  temperature  waveform  to 
have  sharper  crests  and  troughs  than  the  ocean  waves. 

E.   PERIOD  3  RESULTS 

Hydrostatic  conditions  throughout  this  period  were  stable 
with  a  Z/L  value  of  approximately  -0.12.   Results  for  this 
period  were  quite  similar  to  those  presented  for  subsets  2 
and  3  of  period  1,  Figures  17  and  18.   The  latter  subsets 
also  occurred  during  more  stable  conditions. 

Since  all  three  subsets  of  this  period  were  similar, 
only  subset  1  will  be  discussed  in  detail.   Results  for  sub- 
sets 2  and  3  are  in  Appendix  A  as  Figures  36,37,38   and  39. 

The  JDF  of  u  and  T  and  the  CMF  of  wave  height  appear  in 
Figure  28A.   The  CMF  zero  isoline  in  Figure  28A  is  rotated 
approximately  70   CCW  from  potential  flow  results,  Figure  15. 
Phase  amplitude  results  in  Figure  29  indicate  the  above 
rotation  is  due  to  the  temperature  waveform  being  shifted 
approximately  45   from  the  wave's  trough.   This  is  to  be 
expected  in  view  of  previous  results  when  thermal  stability, 
in  terms  of  Z/L  values,  was  considered.   During  period  1, 
Z/L  values  varied  from  -0.03  to  -0.20  to  -0.08  and  maxima  in 
the  temperature  phase  amplitude  varied  from  180   out  of  phase 
with  the  wave  trough,  Figure  23,  to  approximately  in  ph  ;e 
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Figure  24.   JDF/CMF  Results  for  period  2  subset  1 
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Figure  25.   JDF/CMF  Results  for  period  2  subset  2. 
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(A)  (B)  (C) 

Figure  26.   Phase  amplitude  results  for  period  2  subset  1 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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Figure  27.   Phase  amplitude  results  for  period  2 
subset  2  (a)  small  waves,  (b)  large  waves,  and  (c) 
all  waves . 
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with  the  wave  trough,  Figure  20.   Period  3, however,  has  a 
consistent  Z/L  value  of  about  -0.12.   Therefore,  phase  dif- 
ferences between  temperature  and  waves  during  period  3  were 
between  the  extremes  noted  in  period  1.   This  can  be  readily 
observed  by  a  comparison  of  the  phase  amplitude  results  for 
T  and  waves  in  Figure  20  and  23  to  those  in  Figure  29. 

With  respect  to  wave  amplitude  classes,  small  and  large 
waves  appear  to  have  the  same  deviation  from  potential  flow. 
This  is  most  evident  by  a  comparison  of  the  wave  height  CMF 
zero  isoline  in  Figure  28A  and  in  Figure  15.   Recall  that 
during  the  changing  stability  associated  with  period  1,  only 
the  small  wave  amplitude  classes  had  a  tendency  toward 
potential  flow,  Figure  18A.   This  is  in  contrast  to  the 
results  during  the  near  neutral  conditions  during  period  2 
where  small  waves   had  a  tendency  toward  potential  flow, 
Figure  24A. 

Panel  28B  and  28C  results  show  the  CMF  variables  T  and  w 
respectively  compared  to  the  JDF  of  the  wave  height  and  its 
derivative.   CMF  zero  isolines  in  both  panels  are  shifted 
away  from  potential  flow,  Figure  15.   In  addition,  panels 
28B  and  C  display  fewer  small  perturbation  than  similar  re- 
sults obtained  in  period  1,  Panel  B  or  C  of  Figures  16,  17, 
18   arid  19  or  period  2,  Panel  B  and  C  of  Figures  24  and  25. 
This  again  was  probably  associated  with  the  turbulence  occur- 
ring during  the  period. 
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Figure  28.   JDF/CMF  Results  for  period  3  subset  1  (a) 
small  waves,  (I))  large  waves,  and  (c)  all  waves. 
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Figure  29.   Phase  amplitude  results  for  period  3  subset  1 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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P.   PERIOD  .4  RESULTS 

Results  are  similar  to  those  in  Period  3.   The  value  of 
the  thermal  stability  parameter,  Z/L,  averaged  -0.12  in  this 
period  but  approaches  neutral  conditions  near  the  end  of  the 
period . 

The  period  ended  with  a  Z/L  value  of  -0.05.   Even  though 
the  stability  decreased  during  period  4  the  results  were 
quite  similar  to  those  obtained  for  period-  3.   However,  the 
near  neutral  conditions  did  not  occur  until  the  very  end  of 
the  period   and,  thereby,  suggest  that  some  time  lag  occurred 
after  the  long  interval  of  stability  in  period  3  before  the 
wind  wave  coupling  adjusted  to  the  new  environment  and  the 
analytical  method  used  could  detect  charges.   Due  to  the 
similarity  of  results  for  this  period,  only  the  first  of  the 
3  subsets  will  be  used  for  discussion.   Results  for  the  last 
two  subsets  are  presented  in  Appendix  A,  Figure  40,  41,  42 

and  43. 

Results  for  this  period,  shown  in  Figure  30,  are  similar 
to  those  for  period  3.   For  example,  in  Figure  30A,  the 
temperature  maxima  is  45   from  the  wave's  trough  as  was  the 
case  during  period  3,  Figure  28A.   Phase  amplitude  results 
are  shown  in  Figure  31. 
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Figure  30.   JDF/CMF  Results  for  period  A  subset  1 
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Figure  31.   Phase  amplitude  results  for  period  4  subset  1 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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VII.   CONCLUSIONS 

This  study  has  been  an  investigation  of  turbulence  data 
obtained  during  near  neutral  to  stable  conditions.   The  tem- 
perature fluctuations,  were  described  with  respect  to  the 
wave's  influence  over  natural  ocean  waves.   The  statistical 
procedure,  JDF/CMF,  was  used  which  made  possible  interpreta- 
tion of  nonlinear  properties  of  wave  induced  temperature 
fluctuation. 

The  general  characteristic  of  observed  temperature  wave- 
forms leads  to  the  following  conclusions. 

1.  The  most  important  property  shown  is  that  waves 
influence  a  scalar  field  such  as  temperature  as  high 
as  8  meters.   Therefore,  observational  studies  con- 
ducted over  ocean  waves  must  consider  the  wave's 
influence  on  a  scalar  field. 

2.  The  suggestion  by  Volkov  (1969)  that  a  temperature 
observation  at  a  fixed  level  above  a  wave  surface  is 
composed  of  a  mean,  a  fairly  well  organized  fluctua- 
tion due  to  waves,  and  a  fluctuation  due  to  turbu- 
lence appears  to  be  valid.   All  periods  show  small 
fluctuations  superimposed  on  a  more  uniform  sinusoidal 
wave  induced  fluctuation. 

3.  The  phase  shift, of  temperature  away  from  wave  crests 
is  0-135   with  an  average  about  84   for  all  periods. 
This  is  in  contrast  to  the  180   phase  shift  prr  i.  cted 
by  potential  flow. 
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4.  The  temperature  profile  leads  the  vertical  velocity 
although  it  may  have  a  maximum  of  45   on  either  side 
of  the  wave  crest.   This  is  in  contrast  to  the  90 
shift  predicted  by  potential  flow. 

5.  As  stated  in  the  background,  Stewart  (1969)  did  not 
consider  waves  capable  of  influencing  a  scalar  tem- 
perature field.   Therefore,  he  considered  vertical 
velocity  and  T  phase  angles  near  zero  or  180   indica- 
tive of  mostly  turbulence  and  angles  +  90   indicative 
of  waves.   This  study  has  shown  that  under  stable 
conditions  neither  extreme  is  attained.   The  reason 
for  not  consistently  being  at  either  extreme  is 
possibly  due  to  the  wave  induced  fluctuation. 

6.  It  also  appears  for  this  one  study  that  if  the 
atmosphere  is  neutral  or  very  stable  for  an  extended 
period,  sea-air  coupling  will  tend  toward,  but 
probably  never  reach,  potential  flow  characteristic. 
In  between  these  two  extremes,  the  greatest  deviation 
from  potential  flow  will  be  seen. 

This  study  shows  that  JDF/CMF  statistical  technique  is 
very  promising  and  should  be  applied  to  other  data  obtaii ed 
during  the  same  period.   Such  a  study  would  determine  if  the 
conclusions  of  this  study  are  consistent.   Unstable  periods 
should  be  included  if  available. 

One  possible  application  presented  for  this  type  of 
study  was  its  use  in  predicting  laser  propagation.   With 
more  studies  of  this  type  it  should  become  possa      ' 
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predict  the  temperature  fluctuations  under  different  thermal 
stratifications.   Therefore,  with  a  knowledge  of  the  thermal 
stratification  and  the  wave  field  below,  a  more  accurate 
refractive  index  corrected  for  the  wave's  influence  could  be 
provided  to  laser  users. 
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APPENDIX  A 
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Figure     32.       JDF/CMF    Results     for    period     2     subset     3. 
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Figure  33.   J D F / CMF  Results  for  period  2  subset  A 
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Figure  34.   Phase  amplitude  results  for  period  2 
subset  3  (a)  small  waves,  (b)  large  waves,  and 
(c)  all  waves. 
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Figure  35.   Phase  amplitude  results  for  period  2 
gubset  4  (a)  small  waves,  (b)  large  waves,  and 
( c )  all  waves . 
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Figure  36.   JDF/CMF  Results  for  period  3  subset  2. 
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Figure  37.   JDF/CMF  Results  for  period  3  subset  3 
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Figure  38.   Phase  amplitude  results  for  period  3  subset  2 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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Figure  39.   Phase  amplitude  results  for  period  3  subset  3 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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Figure  40.   JDF/CMF  Results  for  period  4  subset  2. 
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Figure  41.   JDF/CMF  Results  for  period  4  subset  3. 
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Figure  42.   Phase  amplitude  results  for  period  4  subset  2 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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Figure  43.   Phase  amplitude  results  for  period  4  subset  3 
(a)  small  waves,  (b)  large  waves,  and  (c)  all  waves. 
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